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Abstract 

A general formlulation for discrete-time quantum mechanics, based on Feynman's method in ordinary quan- 
tum mechanics, is presented. It is shown that the ambiguities present in ordinary quantum mechanics (due 
to noncommutativity of the operators), are no longer present here. Then the criteria for the unitarity of the 
evolution operator is examined. It is shown that the unitarity of the evolution operator puts restrictions on 
the form of the action, and also implies the existence of a solution for the classical initial-value problem. 
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O Introduction 

Discrete-space-time physics is an old tradition originated in solid-state physics. This has been the starting 
point for lattice physics. On the other hand, continuous-spce-time field theory has still some problems; not 
only from the mathematical point of view, but also from the computational point of view. The first kind 
of problems is, essentially, the lack of an exact definition for functional integration, and the problem of 
ultraviolet divergences in interacting field theories ([1], for example). The second kind of problems is related 
to the fact that most of the numerical results of field theory, are in fact perturbative results. A promising 
answer to this kind of problems seem to be lattice field thories, specially lattice gauge theories [2]. These 
theories, however, are interesting, not only as an approximation for continuous-space-time theories, but also 
as independent models [3-7]. 

Another problem, very much related to the above problems, is the lack of a consistent theory of quantum 
gravity. In the context of quantum gravity, there arises a natural scale for space-time, the plank scale, and it 
seems plausible that something new should happen at this scale. In string theories, this scale is the size, or 
the tension, of the string [8]. One can also use this scale as the size of a possible space-time lattice, in fact, 
there is no true reason for the continuousness of space-time: all of the measurements of space-time (direct 
or indirect) have a certain resulotion, which is very much larger than the plank scale. On the other hand, 
there are examples for theories which force the time to be discrete ([9,10], for example). 

Also, there has been attempts to discretize the time in quantum mechanics ([11], for example). In 
this paper, first a general formulation of discrete- time quantum mechanics is introduced (section I). This 
formulation is based on the action principle, or the Feynman's path-integral formalism. 

In section II the unitarity of the evolution operator is exploited to deduce restrictions on the form of 
the action. This section addresses a very old problem of classical mechanics: the problem of nonequivalent 
Lagrangians which give rise to same equations of motion [12,13]. It is shown that in one-dimensional space, 
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all possible forms of the Lagrangian are equivalent to the difference of a kinetic term and a potential term, 
where the potential is an arbitrary function of the position. This means that, at least in one dimension, 
any two Lagrangians (fulfilling the unitarity criteria and) giving the same equation of motion, must be 
equivalent, up to a constant multiplicative factor. At the end of this section, an example of an allowed action 
is presented, which gives rise to the Lagrangian of a charged particle moving in an electromagnetic field. 

In section III, the equation of motion of the Hesinberg position operator is investigated. It is shown that 
this equation is, unambiguously, deduced from the classical equation of motion. In continuous-time quantum 
mechanics, this is not, generally, the case; that is, one can not , in a straight forward manner, deduce the 
equation of motion of the operators from the classical equation of motion. At the end of this section, it 
is shown that the unitarity criteria also guarantee the existence of a solution for the classical initial-value 
problem. 

I Formulation of discrete-time quantum mechanics (Schrodinger picture) 

Ordinary (continuous-time) quantum mechanics is based on two general parts: kinematics, which involves 
the definition of observables, state space, and measurements; and dynamics, which discusses the evolution of 
states. We need not change the kinematics. As for dynamics, we accept that a quantum system is a linear 
dynamical system, the evolution of which is a unitary one. In ordinary quantum mechanics, these among 
whith the principle of correspondence lead to the Schrodinger equation: 



In the case of discrete time, however, we do not have an infinitesimal time-translation. So the Hamiltonian, 
which is the generator of time-translation, does not arise naturally. One can of course write the single-step 



(1.1) 



evolution operator U as 




(1.2) 
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But if -ff is a well-defined Hermitian operator, this is not a real discretization of time; this is only a discrete 
sampling, because one can also define U {t) for arbitrary t as 



JJ{t) = exp 



itH 



(1.3) 



and use it to find the solution of the Schrodinger equation (I.l). 

There is an alternative which does not suffer from these artifacts, and that is to use the Feynman path- 
integral formulation. One only needs to translate the concept of path to discrete case, which is obvious. In 
ordinary quantum mechanics we have 



U{x',x"]t',t")=A j [Dx (t)]exp|^5[a;(i)]|; 



(1.4) 



where S is the action, 



U{x',x";t',t") := {x'\U{t',t")\x"), 



(1.5) 



and I a;' > is the eigenvector of the position operator. Integration is over all paths which satiafy the boundary 
conditions 



x{t') = x', and x{t") = x" 



(1.6) 



This means that we are dealing with a "multiple" integral, whose measure is (rather heuristically) 



[Dx{t)]= H dx{t). (1.7) 
t'>t>t" 

Although this description is rather ambiguous in continuous time, and must somehow be regularized, it 



is completely clear in discrete time. In this case we have 

Un\n"{x',x") = An',n" j dXnexp 

n'>n>n" 

Equivalently, one can use the single-step evolution operator 



n 



S {p^n' ) * * * ; ^n" ) 



(Li 



Un+l/2{x\x") := Un+l,n{x' ,x") = ^„+i/2 Cxp 



n 



S{x\x") 



(1.9) 



We take this to be the axiom of time evolution, substituting (1.2). In (1.8) and (1.9), A is a constant 
independent of x' and x" . 

Now, as we are deahng with a dynamical system, we have 

Un' ,n'"{x' ,X"') = J dx"Un',n"{x',x")Un",n"'{x",x"'), (1.10) 

which leads to a familiar result for the action: 

Sn',n" [x] + Sn",n"' N = Sn' ,n"' [x\- (I.ll) 

To summerize, we define an action for a unit-time interval as 

Sn+i/2{xn+i-,Xn) := 5„+i,„(a;„+i, a;„). (1.12) 

Then, we generalize this to arbitrary intervals as 

Sn',n"{Xn',--- ,Xn") := ^ Sn+l/2{Xn+l,Xn)- (1.13) 

The single-step evolution operator is then 

U{x, y) := Aexp 



n'>n>n" 



{x,y) 



(1.14) 



We have droped the explicit time-dependence of the action for simplicity. But, as one can readily see, this 
does not alter our later results. 

Until now, we have not exploited the unitarity condition of U, and we have no restriction on the action. 
In the following section we will use this condition to restrict the form of the action. 
II Unitarity and restrictions on the action 
The unitarity condition for U is 

J dyU{x, y)U* {z, y) = S{x - z), (II.l) 



or, in terms of the action, 



dyexp <{ - [S {x,y) - S {z,y)] \ = 5{x - z). 



(11.2) 



Expanding the exponent as 



S{x, y) - S{z, y) = {x-z)- V,5(x, y) + 0{\x - z\^), 



(II.3) 



and defining 



u := -{x-z), 



(11.4) 



it is easy to see that (11.2), to lowest order in h, gives rise to 



\A\^ j dyexp 



-{x-z}- \/^S {x,y) 



5{x — z). 



(11.5) 



The action, itself, may depend on h. So, to be more exact, the above equation holds for zeroth order (of h) 
term of the action, i.e. the classical action. But, as we are going to restrict the form of the classical action, 
this does not mind. 
Defining 



Y~ V,S{x,y), 



(11.6) 



we have 



dy 



det 



dy 

dY 



exp 



n 



{x~z)-Y 



5{x — z). 



(11.7) 



Now, the Jacobian is a function of Y and x only: It does not depend on x — z. So the left-hand side of 
(II. 7) is the Fourier transform of this Jacobian, and is supposed to be the Dirac delta distribution. The 
Jacobian should, therefore, be independent of Y. Since A is independent of x, the Jacobian should also be 
independent of x. We therefore have 



det 



dy 
dY 



const.. 



(II.) 



det ( ^ ) = const. (II.8)' 



or 

This is a nontrivial differential equation, which restricts the form of S. If the dimension of the space is 
more than one, it is not easy to determine the most general solution of this equation. But in one dimension 
the life is simpler and one can obtain such a solution. 

II. I Action in one-dimensional space 

In one dimension, we have to solve 

- (II.9) 



dxdy T ' 

where we have written the constant as above for later convenience, r is the time step parameter. This 
equation can easily be solved: 

777 

S{x,y) = xy + fix)+g{y), (11.10) 

T 

or 

TTi T 

S{x, y) = -{x-y)'-- [V {x) + V {y)] + {x) - ^ {y)] , (11.10)' 

The third part of this expression has no effect on the dynamics of the system. It is easy to see that addition 
of such a term to the action is equivalent to the following gauge transformation 



\x >— »■ exp 



\x > . (11.11) 



It is also easy to see that this term does not affect the equation of motion (which we will encounter in section 
III). We therefore conclude that the most general solution of (II. 9) is equivalent to 

S{x, y) = ^(x-yr-^ [V (x) + V {y)] (11.12) 

This form is a very familiar one. In fact, if we divide the action by r and let r tend to zero, we obtain 

L [X it)] := lim ^ = fx^ - V{x) (11.13) 

r— >0 T Z 
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which is the standard form of one-dimensional Lagrangians used in textbooks. This result is, however, an 
important one; First, in classical mechanics the general Lagrangian formulation does not restrict the form of 
Lagrangian as a function of position and velocity; one can have terms in Lagrangian which are, for example, 
quartic in x. However, nobody has ever needed to consider such Lagrangians for real one- dimensional systems 
and it seems that nature has chosen the special form (ILIS). This formulation provides an explanation for 
this fact, based one the general assumption of unitarity. Second, there are many nonequivalent Lagrangians 
that lead to a same equation of motion [12,13]. (By nonequivalent Lagrangians we mean Lagrangians, the 
difference of them is not a total derivative.) One essentially has no way to choose one of this Lagrangians 
for quantum theory. The discrete-time formulation characterizes only one single Lagrangian among these. 
II. II Action in multi-dimensional space 

It is easily seen that an action like (n.l2) satiafies (ILS)'. To have a taste of other kinds of solutions, 
consider a perturbative approach: 



S{x, y) = So{x, y) + s{x, y) := ^{x -yf-\ [V {x) + V {y)] + s{x, y). 



(11.14) 



We want to solve the linearized equation, corresponding to (IL8)', in terms of s. This equation is readily 

seen to be 



\ dxdy J 



(11.15) 



A special solution of this equation is of the form 



S= 2{[-^l (a;i,y2,2/3,---)--^l {yi, X2, X3, ■■■)] + [-^^2 {yi, X2, X3, ■■■) -A2{Xi,y2, X3, ■ 



+ •••}. (n.i6) 



In the limit r ^ 0, we have 



s = 



dAi . 
dxi 



dAi . dAi . 

X2 + — X3 + ■ ■ 



dX2 



dxr. 



dA2 . 

^ ^2 

OX2 



dA2 . dA2 . 

X\ -h — X3 -h ■ 



dx\ 



3x3 



, (11.17) 



or 



'dAi . dA2 . 

T I — Xl + — X2 + ■ 



dxi 



dXn 



T / dAi dA2 
2 \~dr'^ "dT 



(11.18) 



The second term is a total derivative. Defining 



we see that this solution corresponds to the Lagrangian 

which is the Lagrangian of a charged particle moving in the potential V{x) and a magnetic field, the vector 
potential of which is A{x). 

Ill Heisenberg picture, relation to classical mechanics, and existence theorem 
for the solution of the equation of motion 

Evolution of the operators in Heisenberg picture is just like the case of continuous time. We first define 
Heisenberg operators at time n in terms of their Schrodinger counterparts as 

:= {U^rh'^U^ (HI.l) 

Now, we want to discuss the equation of motion for these operators; to be more specific, we want to show 
that the Heisenberg position operator satisfies the "classical" equation of motion. In fact, the true classical 
equation of motion is (or can be) somehow ambiguous for operators, because it involves the position operator 

at three distinct times, and these operators do not commute with each other. But we will see that a certain 

ordering is dictated from the evolution governed by (HI.l). 

„ . , , . , ,1 dS(x\w) , ,, , , „ 

Consider the matrix element < x \ \x >, where we define 

ow 

:= U^xU, and x^ := UxU'' . {111.2) 



We have then 



< x'\^^^\x" > =< x'\U^^^^U\x" > 
ow ow 



= Jdz< x'\U^^^^^\z > U{z,x") . (ni.3) 
= J dzU*{z,x')^^^^U{z,x") 



Now, setting w = x", we come to 

where the momentum operator p is the generator of space translation, and its matrix element satisfies 
equation (III. 4). To deduce (III-4) from (III. 3), we have exploited the form of U in terms of the action, and 
also unitarity of U. 

Equation (III. 4) can be rewritten in terms of operators themselves: 

,,dS(x\x") , II i,-^dS{x^ ,x) II 

^ dS{x\x) _ . 

o o = -p, (III.5) 

where the left-hand-side of (III. 5) is an ordered form, in which is are at the right of x^'s. 

By a similar arguement, or by taking the Hermitian conjugate of (III. 5), one concludes that 

^^^^ = -p, (III.6) 

where the left-hand-side of (III. 6) is in opposite ordering. Finally, since in this case both of these orderings 
lead to the same result, we can define a nonoriented ordering, which is equal to both of the above orderings: 

dS{x\x) 



o- 



dx 

By a similar arguement we come to 



-o = -p. (III.7) 



dx 

Now, eliminating p from (III. 7) and (III. 8), we obtain 

dS{x,x^) dS{x\x) 



dS{x,x^) , 

o — o=p. (III., 



o- 



O + O ^-^0 = 0, (III.9) 



dx dx 



or 
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But this is the classical equation of motion for position operators (obtained by extermizing the action), 
except that it is time ordered (either in the forward, or in the reversed direction). It has no ambiguity and, 
once we know the classical action, everything is determined. 

Now we come to the problem of existence of solution in classical mechanics. Suppose that we begin by 
an arbitrary action, which does not satisfy the unitarity criterion. One can not guarantee that the classical 
equation of motion has a solution for Xn+i for any choice of a;„ and Xn-i- this equation may be a nonlinear 
complicated one. This means that, not every initial- value lead to a path: there may be a time when the 
particle can not go anywhere. 

Now suppose that the unitarity condition holds, so that the equation of motion for the operators, (111.7) 
through (III. 9), hold. For any pair of initial values xq and X-i, one can solve the classical counterpart of 
(III. 8) for po. Use a Gaussian wave packet with 

<x>=xo, <P>=Po, (III.ll) 

and 

Ax = a^, Ap = a-'^, (III.12) 

as an initial state for quantum mechanics. Then consider equation (III. 9) in the limit h ^ 0. In this limit, 
the uncertainities tend to zero, the operators commute, and we have 

dS{< x>,<x^ >) dS{< x^ >,<x>) ^ 

d<x> d<x> ■ y ■ ) 

Notice that there always exists a "< >", because one can compute it through the Schrodinger picture. 

Now, we know that 

os(^ = ,3 =< ^ >= ^^«;>:<^^», (111.15) 

OXo o < x > 

where the last equality holds only in the limit % ^ 0. Substituting xq and X-\ iov < x > and < x^ > in 
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(III. 16), we will have 

OXq OXq 

which means that < > (in the limit ?i — > 0) satisfies the classical equation of motion for the initial-values 
xq and x-i. Specially, if the dimension of space is one, we have seen that any action which satisfies the 
unitarity condition is of the form (11.10)', or equivalently (11.12). The equation of motion for such an action 
is linear in a;„+i and the coefficient of a;„+i is nonozero. Therefore every initial-value problem has one and 
only one solution. 

In general, we have proved that the unitarity of the evolution operator, constructed through (1.14) from 
the action, is sufficient for the existence of a solution for classical initial- value problems. 
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